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SEMI-ANNUAL STATUS REPORT - NSG 5014 - JUNE 30, 1977, Dartmouth College 
The D^u:tlIlouth LANDSAT Research Group, working in collaboration with tiie Goddard 
Institute for Space Studies staff, continued the study cuid analysis of LANDSAT data 
together with appropiate field infonaation (ground truth). Technical development 
continued to improve use of the avail^dDle data processing devices, especially the 
RAMTEK display. A user's manual was completed. 

1. Ag ricultural Crop Studies (Prof. Huke) 

Without graduate degree candidates or full-time research assistants, the 
principle technical efforts in this area of interest have been ad hoc undergraduate 
student projects which result in limited, in-house reports or execution of specific 
tasks undertaken at the request of GISS staff. 

a. Some effort has been applied to the data available on the contents of the entire 

9000 acre TWIST Plantation, but no conclusive results are available. However, aerial 
photos, 512 chcmnel spectrometer data, LANDSAT data and much ground truth 
are on file for analysis. 

b. Considerable effort by several different students was put l\ito the so-called 
Illinois (Peoria) Project, in an attempt to classify corn and soy crops from 
LANDSAT data. Although a number of different approaches were made to this 
problem and many aspects of the solution were examined, it remains an open- 
ended project. 

c. A number of short student projects have been accomplished in the field of 
agro-meteorology with Professor Malmstrom. 

2. Forestry (Prof. Huke) 

Endly Bryant is the principal researcher working in association with Gibb Dodge, 

RC & D Area Forester with the Cooperative Extension Service, University of New 
Hampshire, based in Conway, NH and with and undergraduate assistcmt. The goal is 
to use conputer classification of LANDSAT data to make a product useful to the field 
forester. 

The main thrust of our work this half year has been a cooperative project with 
Seven Islands Land Company in Bcuigor, Maine . We developed classifications of thrw-e 
Maine townships, and compared the acreages of softwood, mixed vpood, and hardwood 
generated with those from a standard inventory that Seven Islands is having done. 

Ttie first classification, using "rational signatures” yielded acreages different 
enough from the inventory that the results would not be useful to Seven Islands. 

We then revised the classification so there was agreement in acreages in a test 
township, and applied it to two adjacent townships. There was close agreement in 
one township and quite a difference in the other. We believe the difference is due 

to the large amount of timber harvest in the township. We are now in the process of 
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extending the classification over 19 additional townships ^ about 500,000 acres 
total - to see if the discrepancies balance out over a larger area* We have had 
cooperation, constructive criticism and enthusiasm from Seven Islands on this 
project* An idea which seems strong to both them and us is an overall forest 
information scheme integrating standard inventory techniques such as aerial 
surveys with potential techniques such as LANDSAT classifications and computer 
displays of information* 

We have also been working with classification of clearcut areas - locating 
training, test, and applications areas, and gathering ground truth for them through 
overflights of the areas and contacts with the companies doing the cuttings. Aerial 
photos of the areas will be flown in July* 

Bryant and Dodge attended the sy rapes ium on remote sensing sponsored by ERIH (Ann 
Arbor), and found it fruitful in ideas and contacts with other people in the field. 
One concrete idea we gained there and put into action here was making color 
printouts of mapped LANDSAT data by using colored ribbons on our printer. 

As part of a computer course here at Dartmouth, DTSS was programmed to classify 
sample LANDSAT data using four algorithms. 

Preliminary work on spruce budworm detection showed that we need more ground 
truth to proceed. 

While doing classifications, same area two different passes - we found a 
discrepancy of about 1.3% in pixel location between the passes. This is enough 
to foul us up somewhat. 

We have been using the ^'polygon progrcim" which enables us to have the 
computer count pixels v;ithin the Maine tomship boundaries - this is a great 
improvement over previous techniques. 

We gave two remote sensing lecttires at UNH and two at Dartmouth, and made ^ 
display for a workshop for foresters in Hillsboro, NH . A paper has been ac -epted for 
a workshop sponsored by the Society of American Foresters (abstract attached - 
Enclosure 1) , to be hold in Tuscon, Arizona, next January and one was published in 
the Mt, Washington Observatory News Bulletin (also attached Enclosure 2). 

3. Land Use (Prof. Huke and Advisor, Prof. Lindgren) 

The principle effort in this area is by undergraduates doing limited paxt-'time 
projects* These projects have provided the basis for an analysis of a combination 
of LANDSAT data and ground truth (aerial photos) for the New Haven, Connecticut 
urban area which has been a subject of study in the Geography Department for some 
time* These results were scheduled to be presented at the Symposium on Machine 


3 


Processing of Remotely Sensed Data at Purdue in June but an unfortunate 
accident to Professor Lindgren prevented completion of the presentation and 
his attending the conference. However, this material has been organized as 
background for possible future student projects, particularly for the application 
of LANDSAT-C data in the coming year. 

We have also cooperated with Dr. Paul E. Bruns of the University of 
New Hampshire and associates who has analyzed changing land use patterns in 
the Dover, New Hampshire area, using LANDSAT data and aerial photography. This 
work was reported in the local press (Valley News Story, 17 March 1977 - 
attached - see Enclosure 3) and will be the subject of a technical report at a 
later date. 

4. Water Quality (Profs. Reynolds and Drake, Graduate Research Assistant McGirr) 

Spectra has been obtained from Lake Powell, and Green and Round Lakes, NY. 

The latter are known calcite precipitating lakes, emd their spectra show interesting 
similarities with those from Lake Powell. VJe are in the process of analyzing 
these spectra, and those from Lake Tahoe and Lake Champlain, in order to develop 
a qualitative discriminator for calcite precipitation in hard water lakes. 

The analysis will be made using spectra stored in the Dartmouth Time Sharing 
System. McGirr spent the winter term retrieving these spectra from GISS and 
transposing them to the Dartmouth system. The data transfer was completed at the 
end of the spring term (May 1977) . Vfe anticipate early completion of the analysis 
and publication of the results soon after%i;ard, unless the findings prove trivial. 

The earlier work by Carolyn Merry on "The Correlation and Qucuitification 
of Airborne Spectrometer Data to Ground Turbidity Measurements at Lake Powell, 

Utah” (Masters Thesis Sept 1976) has now been published as US Army Corps of 
Engineers CRREL Report, October 1976. 

5. Mineral Resources (Prof, Birnie, Research Asst. Hutton and several undergraduate 
assistcints) 

A. Reflectance SPectra of Vegetation Growing on Mineralized Ground 

Data presented at the Fall 1976 G.S.A. meetings by Birnie and Hutton showed 
that spectra collected and averaged within a mineralized zone differed from those 
collected and averaged outside of the mineralized zone. The spectral difference 
is concentrated in the 0.67 region of the reflectance spectra and is clearly 
manifested in a channel by channel ratio plot of "mineralized/nonmineralized" 
spectra. Collection of tliese averaged spectra required prior knowledge of where 
the mineralized ground is located. 
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Since December, 1976 progress has been made in discriminating mineralized 
vs nonmineralized ground without prior knowledge of the geology. The procedure 
involves 4 steps: 

1) isolation of uniformly dense tree spectra with a rejection of 
spectra contaminated by grass, roads, etc. This is accomplished 
by a discriminant in the 0.67 ^ region of the spectra. 

2) separation of dense tree spectra into shaded and sunlite categories 
This is accomplished by a discrimincint in the 0.79 yU region of 

the spectra. 

3) separation of mineralized and nonmineralized spectra for each 
of the groups of shaded and sunlite spectra. This is 
accomplished by a discriminant in the 0.67yu region of the spectra. 

4) recombining of shaded and sunlite spectra and making a map showing 
regions of mineralized and nonmineralized spectra. 

This work which was done using spectra obtained with multi- channel 
spectrometer flo\7n in an aircraft, provides a basis for recommendations 
concerning preferred spectral characteristics of satellite sensors and possible 
methods of applying existing and future LANDSAT data to mineral exploration. 

A paper on this work is in preparation and early publication is anticipated. 
B- Reflectance Spectra of Exposed Hydro the3nnal Altered Rock surfaces using 
LANDSAT Data (Prof. Birnie, Graduate Research Asst- Dykstra) 

A computer analysis of LANDSAT data from the Chagai Hills area of 
Pakistan has been undertaken - 

A field visit was made by Jon Dykstra to the Chagai Hills on October- 
December, 1976 to obtain preliminary geologic ground truth necessary for 
the computer study. 

As part of the computer study, Dykstra has developed a modified 
classification algorithm, based on principal component analysis, for 
processing the LANDSAT data applicable to the detection of possible copper- 
bearing areas. This part of the work will be reported at the Geological 
Society of America Annual Meeting in Seattle, November 7-9, 1977 (cibstract 
attached - Enclosure #4). 

(A further visit to Pakistan is planned by Birnie and Dykstra - 
to be funded by the National Science Foundation - during the late 
summer 1977, to check out on the ground various anamolies which have shown 
up in the geologic maps in the classification work.) 
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C. LANDS AT Satellite Mapping of Ultramafic Rocks (Prof. Stoiber cind 

Graduate Research Asst* Carlson) 

Work involving the detection of ultramafic rocks -- usually associated 
with copper ore ^ using LANDSAT data, in Guatemala and Costa Rica was 
reported at the 5th Reunion of Central American Geologists in Managua, 

Nicaragua, on 24 February, 1977 (Abstract attached ~ Enclosure #5) . 

Work continued in mapping the Oman ultreimafics cind in refining a statistical 
signature definition and mapping procedure- The existence of gossan over 
massive pyritic copper deposits was learned from a mining company operating in 
the area, and definition of these gossans with the satellite data will be 
attempted- Separation of ultrcimafic rocks from others in the area has been very 
successful and in connection with these studies, Carlson has developed an 
improved algorithm for application of LANDSAT data analysis to these specific 
geologic problems and perhaps in other fields- (This improvement is being 
reported at the Pecora Symposium, Sioux Falls (ND) , October 30 - November 2, 

1977 (abstract attached “ Enclosure #6) and at the GSA Annual Meeting, 

Seattle, November 7-9, 1977 (abstract attached - Enclosure #7). 

6. Remote Sensing of Volcanic Emissions (Prof. Stoiber, Research Asst. 

Bratton and Graduate Research Asst. Malinconico) 

The Dartmouth research group v/orking on remote sensing of volcanic 
emissions continued its analysis of data obtained earlier from volcanoes 
in Central America, Italy, Sicily, Hawaii and continued modification and 
timing of instruments (SO^ and HCl sensors) to make them more field-worthy. 

Work concerned with the Central American volcanoes was reported at the 
5th Reunion of Central American Geologists, Managua, Nicaragua, 24 February 
1977 (abstract attached - Enclosure j}8) . 

Part of this work was also incorporated in the Masters Thesis of 
Richard W. Hazlett, Dartmouth College, June 1977 "Geology and Hazard of the 
San Cristobal Volcanic Complex, Nicaragua" which utilized remotely detected 
data (SO^) by the project COSPEC. In addition, Malinconico is conducting 
a study of lava flow patterns as depicted by LANDSAT data to supplement and 
improve volcanic maps developed by Hazlett. 

In March 1977, Stoiber, Bratton and Malinconico undertook further field 
tests of the COSPEC and GASPEC in vicinity Mt. Etna, Sicily, An important 
discovery is the problem of dependence on the geometry of sight lines for 
the reception of SO^ spectral data. HCl (GASPEC) sensing has not been 
satisfactory and the machine is still being analyzed. (Trip Report attached - 
(Enclosure #9) . 


Bratton field tested the COSPEC in Alaska in April 1977 but results were 
inconclusive and warrant further comparison with airborne detection data obtained 
by a University of Washington party • 

Progress is underway on developing a smalls remote controlled aircraft sampler 
for providing simultaneous air sampling data in coordination with COSPEC/GASPEC 
tests . 

At the present time, concentrated work on bringing the GASPEC (HCl) instrument 
to acceptable stability and sensitivity is planned as well as a further visit to 
Mt. Etna summer 1977- 


SATELLITES FOR PRACTICAL FOREST MAPPING? 

/ by 

A, Gibb Dodge 

RC & D Area Forester 
and 

Emily Bryant 
Dartmouth College & 

Goddard Institute For Space Studies 

After World War I, forest managers began using aerial photos 
as a supplement to cruising in gathering forest information. Since 
then, aerial coverage has proved itself an essential part of forest 
inventory. Once more, we have a new source of forest information - 
satellite data. We believe that computer classification of satellite 
data can fit into forest inventory systems as aerial photo 
interpretation has. 

Two satellites, LANDSAT 1 (launched in 1972), and LANDSAT 2 
(launched in 1975) record the earth's reflectance acre by acre from 
a height of 550 miles. They transmit the data to earth where it is 
put on computer- compatible tapes. Since different land cover types 
usually reflect differently, we can train a computer to read, the 
tapes and recognize land cover types by their reflection. A computer 
printout of the satellite data gives a 1:24,000 scale map classifying 
each acre into one of up to 30 categories. 

As a test, Seven Island Land Company, Bangor, Maine, agreed to 
compare detailed traditional forest inventory data on three tov/nships 
(57,955 acres - 23, 454 hectares) under their management with acreage 
tallies and maps derived from computer classification of data from 
LANDSAT 1 (scenes 5480-14043 and 5444-14073) . Digital maps including 
three forest categories, water and open area were produced using 
techniques developed earlier in northern New Hampshire by Dartmouth 
College, Goddard Institute for Space Studies, and Cooperative Extension 
Service, University of New Hampshire. Aerial photography and an 
overflight of the three townships in Maine provided us with ground 
truth . 

Preliminary comparison of satellite data with Seven Islands Company 
data indicates that satellite data consistently underestimates softwood 
forest type. Additional analysis and collaboration with Company 
forest managers will make LANDSAT digital maps more useful in 
measuring and locating general forest types for tax purposes in the 
test area. 


(For presentation at Integrated Inventories of Renewable Natural 
Resources Workshop, January 1978, Tucson, Arizona.) 
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i GEOLOGIC APPLICATIONS OF COMPUl'ER PROCESSED LA^IDSAT DIGITAL DATA; 

CttAGAl HILLS, BALUCHISTAIH, PAKISTAN 

li DYKSTlUi, Jon D,. , and BIRNIE, Richard VI. , Department of Earth 

Sciences, Dartmouth College, Hanover, Mew Hampshire 03755 
• A geologic map of the Chagai Hills has been generated by computer 

'} processing of lANDSAT digital data, ^he classification algorithm, 

|| based on principal component analysis, was developed at Dartmouth in 
i collaboration with the Goddard Institute for Space Studies, NY. A 

principal component vector in the LANDSAT 4-diraensional color space is 
calculated from the pixels of control areas representing each of the 
i geologic types to be mapped. An ^mknown pixel is classified as the 

" geologic type whose principal component vector passes closest to it in 

I the color space. A confidence value for the classification of each 

pixel is calculated. This value increases with 1) the pixel’s 
5 proximity to the principal component vector and 2) the separation of 
I tlie different principal component vectors witliin the color space. 

I Spurious pixel variations due to inho.mogenities within a pixel are 

I smoothed by averaging each pixel with its eight contiguous neighbors. 

The algorithm v/as applied to 5000 km^ of the Chagai Hills. This 
I region is located in the southeastern end of the Mediterfanecin-Iran 

! copper belt and is a likely porphyry copper setting. Cretaceous 

f andesite flows and pyroclastics are intruded by Cretaceous to 
1 Eocene (?) plutons ranging in composition from quartr-diorites to 
I granites. The classification procedure distinguished three rock 
\ types (volcanics and two intrusive phases) , four vinconsolidated 

deposits, and zones of possible hydrothermal alteration. LANDSAT 
geologic mapping is especially applicable to the Chagai Hills due to 
[ the region’s large areal extent, arid climate, well exposed geology, 

; and very limited accessibility. 


! 


r 

□ Oral Kl Poster 


□ J3kher 


n Symposiuni ^ 

(tltlj of symposiuiTi) 






CLASS! Pi CATION 
You must specify n/tc. If more 
thJin one category h appropriate, 
ituHcate your order of preference 
hy numbers, He specific. 


* gcochtrmistry 


archeologic 

coal 



education 

engineering 

environmental 

extraterrestrial 

general 

history of 

marine 


Precainbrian 

Quaternary 

structural 


gcornorphology 


geophysics 


geoscience in format Ion 
hydrogeology 

m i n eralogy/ cry *> tal I ograp hy 
pnleomagnetism 

paleontology/paleobolany 


petrology 

e.xperitnental 

igneous 

metamorphic 

sediinentology 

sedimentary petrology 
stratigraphy 


tectonics 


Kernote 


1st 

Sensing. 


Speaker Dy.}:r.tca 0 Student paper □ GSA Student Associate 

I v/illbe available to serve as a cochairman for a technical session on or concerning 

For correspondence purpor^es, list address of senior author if different from above 


Fhorie luiiub'^rs and djics where senior author can be contnetod _ 
Geological Society of America Annual Meeting ^ 



(uO.t) rrospnt L.O Mrxitincf 

Seattle/ Washington November 1-3, 1977 

^ ^ Enclosure ir4 


3/7/ 




'i'o be delivered at the 5th Reunion of Central American. Geologists 


ABSTRACT 

Remote Sensing of Central America Ultramafic 1'ock.s 

Using LAKDSAT Data 

by . 

Gerald G. Carlson . /7/^ 

Department of Earth Science s_ ^ , 'L)cx.r~h^o-,e^/^ 
Hanover/ N.K. 03755/ U.S.A. . > 

A program of remote sensing of natural resources using 
digital IiANDSAT satellite data has been underi'/ay for nearly • 
three years at Dartmouth. One project has involved the ' 

delineation of ultramafic rocks in areas of varying climatic 
conditions. It is well known that in many areas distinctive 
vegetational anomalies are developed over ultramafic rocks due .a 
to the sterility of soils derived from these rocks. 

In Central America, ultramafic rocks occur along the 
Motagua and Polochic valleys in Guatemala and the Santa Elena 
peninsula in Costa Rica. A distinctive lack of vegetation is 
evident over many of these known ultramafic exposures. A 
computer classification of the satellite data over two areas, 
near Huehuetenango in Guatemala and the Santa Elena peninsula, 
has shown that it is possible to distinguish most of the ultra- 
mafic rocks. Areas of false classification, those which are 

mapped with the satellite data as ultramafic .rocks but are not, 

“ .* 

include fresh alluvium in river valleys, vegetated alluvium . 
and dry farm fields in areas of arid climate, and fresh, lightly 
vegetated volcanic flov/s. Further confusion arises with the ■ 

Costa Rica example; burning of grassy areas has occurred both 
on and off the ultramafic rock outcrop just prior to the satellite 
pass . 

5th Reunion of Centrcil Mierican Geologists, Mangua, Nicaragua, 24 Febuary 1977 
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ABSTRACT 


GEOLOGIC t4APPING BY COMPUTER ANALYSIS 
OP DIGITAI. LANDSAT DATA 

CARLSON, Gerald G, , Departraent of 

Earth Sciences, Dartmouth College, 

Hanover, NH 03755 

Digital LANDSAT data is applied to geologic mapping 
at Dartmouth College using an IBM 360/95 computer and 
software assistance at the Goddard Institute for Space 
Studies in New York, together v^ith our own mapping algorithm - 
LANDSAT computer compatible tapes used have been geometrically 
corrected using Bernstein’s method. Maps are printed at an 
accurate 1:24,000 scale on a high speed printer. 

An algorithm has been developed which has tv70 special 
features: assignment of SIGHATURE PROBABILITIES and refinement 

of pixel classification considering NEAREST NEIGHBORS . Four- 
band signatures are extracted using multiple training sites for 
each classification category, each site on the order of 
25 pixels in size. The four-band reflectance values from 
LANDSAT are no ruralized to unit vectors and compared to a 
standard signature set. The probability of each standard 
signature representing each training category is calculated. 
These probabilities are modified by an albedo v;eightlng factor. 

The signatxire reflectance characteristics of a specific 
rock type typically shov? a broad variation . This is due to 
the natural inhomogeneity of the exposed svirface, resulting from 
factors such as v/eathering, soil cover, vegetation cover and 
slope direction and angle. This variation, combined v.'ith the 
noise factor inherent in tlie satellite datci, results in a 
large degree of overlap of spectral characteristics among most 
Pecora Symposium, Sioux Falls (KD) , October 30 -November 2, 1977 
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rock types. To overcome this probleia^ v/e re-calculate the 
probabilities for each pixel v/ithin a survey area by averaging 
it v/ith its eight adjacent neighbors. Root mean square 
averaging is done to emphasize high probabilities. The final 
map therefore delineates areas in V7hich there is a high 
probability of the occurrence of one of the rock types 
originally chosen as a training category. 

The mapping algorithm has been applied to several specific 
geologic problems v^ith good success in delineating the target 
features and significant improvement over earlier supervised 
classification methods. These examples include outlining 
ultramafic rock exposure in Western Guatemala, mapping members 
of the ophiolite sequence in Oman and defining the areal extent 
and relative ages of pre-historic flov;s from Volcan Sem 
Chris tobal in Nicaragua. 




I'LliASE SUHMI’l’ ONE ORIGINAL AND J^OUR COPIES 


j T.!i Geoto^cal Society of America 
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GEOLOGIC MAPPING BY COMPUTER ANALYSIS OF DIGITAL LANDSAT 
DATA 

CARLSON/ Gexald. G. , and. STOIBER, Riciiaird E. / Depairtmanb 
of Earth Sciences, Dartmouth College, Hanover, 

New Hampshire 03755 

Digital LANDSAT data from computer compatible tapes geomet- 
rically corrected by Bernstein ' s method are used for 
geologic mapping using an IBM 360/95 computer and software 
assistance at the Goddard Institute for Space Studies, 
together with our own mapping algorithm. Our algorithm has 
two special features: assignment of Signature Probabilities 

and refinement of pixel classification considering Nearest 
Neighbors. Four-band signatures are extracted using .multi- • 
pie training sites for each classification category.. These 
are normalized to unit vectors and compared to a standard 
signature set. The probability of each standard signature 
representing each training category is calculated and 
modified by an albedo v/eighting factor. Spectral reflec- 
tance characteristics of a specific rock type typically show 
a broad variation due to weathering, soil cover, vegetation 
cover, slope direction and angle and noise in the satellite 
data. This results in overlap of spectral characteristics 
among most rock types. So, for each pixel in a survey area, 
v7e re-calculate its probability of representing each rock 
type by averaging it with .its eight neighbors. Root mean 
square averaging emphasizes high probabilities. The final 
map delineates areas with a high probability of occurrence 
of each rock type chosen as a training category. The 
algorithm, used to map members of the ophiolite sequence in 
Oman, successfully delineates ultramafic rocks, mafic flows, 
older sediments and terrace deposits. 
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GASES AT CENTRAL AI>iERICAl<f VOLCANOES 


by Richard E. Stoiber, Dartmouth College, Hanover, U-S.A. 

Monitoring of temperatures and gases is continuing at 
three volcanoes in Central America, San Cristobal and Cerro Negro 
in Nicaragua, and Santiaguito in Guatemala. The gases are being 
collected as condensates from v;hich relative ’amounts of, HCl, SO 
and HF can be determined. SO^ is also determined by remote sensing 
using a correlation spectrometer. Remote sensing of SO^ has 
indicated a lower rate of emission at San Cristobal since the ' •- 
eruption of March, 1976. Now it is hundreds of tons per day. 

Then it v/as thousands of tons per day. In November, 1971, before 
the present episode, the rate was 360 tons per day, about equal 
to that in November 1976. Santiaguito which erupts ash about once 
per hour in clouds v/hich are five hundred to a thousand meters 
high, produced less SO 2 per day than San Cristobal where no ash 
is being discharged. Measurements in July 1976 indicated less 
than 200 tons per day. 

The remote sensing data seam to be of sufficient value to 
v;arrant consideration of frequent and regular monitoring coupled 
with monitoring of sreismic activity associated with a volcano . 
Regular monitoring must be carried out by a local organization. 


5th Reunion of central American Geologists, Managua, Nicaragua, 24 February 1977 
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. REPORT ON TRIP TO SICILY - March 17-31 

1977 

(NASA Grant NSG 5014) 





Purpose of Trip . To test the remote sensing, HCl instrument. 

To obtain nev/ data regarding sulphur dioxide emission at 
Mt. Etna using the SO 2 spectrometer. Field tests were planned 
. on Mt. Etna where the fumes emanating from the summit can be 
viewed remotely and where the temperature is colder than in 
Central America. It was believed that this would be an optimum 
place for observing HCl. 


Results* . Mt. Etna v/as fuming copiously. We measured the HCl 
V7ith difficulty. The instrument was unstable both electronically 
and mechanically, but it gave us some data. The best data was 
obtained by observing the fumes at the crater and poor data were 
obtained observing fumes some distance from the crater. This 
v?as assumed to be because of lov?er concentrations at greater ' 
distances. It is not certain V7hether these data can be evaluated 
qualitatively, but general answers can presumably he obtained. 

SO 2 measured at several distances from the crater on several 

days was between 1,300 and 1,900 tons/day. The results were 

consistent and thought to be good measures of the SO 2 abundance. 

There v;as no evidence that SO 2 dissolved in v;ater vapor more 

extensively at the slightly greater distances at which some 

measurements were made relative', to others. These results compare 

with the 1975 data of Zettwog, which averaged 3,500 tons/day and 

the earlier data of Stoibar and Malone which averaged 230 tons/day. 

The emission of water vapor appeared visually to be very much greater 

than during the March 1975 visit of Stoiber and Malone. The 

* All rcsullzs aro boing re-cvaluaterl in the laboratory. The data in. this 
report are simply tho.'ie acquired from calculations made in the field. 


The large change factor of probably 8 times v/as not surprising. 
Difficulties in measuring arose from the fact that the plume 
changed direction rather quickly during the day, sometimes during 
a single run. It v/ould seem best to visit in a season in X'/hich 
the v;ind direction was more constant. If indeed meteorological 
data, V7hich vtq have requested, indicates that there is such a 
time of year. 

A camera with motor-drive was field-tested for the first 
time. The time-lapse photography taken with the camera has been 
sent for development and results will be evaluated. 

Future Plans . At present it is planned to check the electronics 
of the HCl machine in detail to improve its stability and possibly 
its sensitivity. It is thought that the machine in the summer 
time could be taken to the summit of Mt. Etna very near the point 
of emission of the plume and v;ould thus probably find very large 
concentrations from which, v/ith stable electronics, we could 
easily calculate tons of HCl emitted per day. The necessity 
for freezing conditions, which had been postulated as a possible 
requirement in using this machine, seems to have proven not 
necessary; a large concentration is the key factor and a stable 
machine is necessary. 

On the same expedition v;hen HCl measurements were made, or 
perhaps another time one could measure the SO^. One could also 
investigate SO 2 emission related to distance from the vent. It 
appears that the configuration of the roads at Etna is such that 
with the constant wind direction and velocity one could make 
several runs per day at different distances from the vent. This 


would allow the evaluation of whether or not there is an 
important amount of the SO^ dissolving in the water vapor 
with time. We nov/ await evaluation of meteorological data to 
determine when v/ould be a most suitable time of year. It 
is contemplated that this may be used as a Master's degree 
project. It is possible that Mr. Malinconico will undertake 
this study. 


Richeird W. Stoiber 



